Distributed source identification for wave equations: 
an observer-based approach 
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Abstract — In this paper, we consider a wave equation on a 
bounded interval wiiere tlie initial conditions are linown (are 
zero) and we are ratlier interested in identifying an unlinown 
source term q{x) thanlis to the measurement output y which 
is the Neumann derivative on one of the boundaries. We use a 
back and forth iterative procedure and construct well-chosen 
observers which allow to retrieve q from y in the minimal 
observation time. 

I. INTRODUCTION 

In a recent work [1], Blum and Auroux proposed a 
new inversion algorithm for identifying the initial state 
of an observable system, based on the application of 
back-and-forth observers. Noting that we have only access 
to a measurement output on a fixed time interval (0,T), 
the idea consists in proposing a first asymptotic observer 
for the system that will be applied in this time interval and 
a second one that will be applied to the system where the 
direction of time is reversed. These observers are then used 
iteratively to get a better estimate of the initial state after 
each back-and-forth iteration. If the two observer gains are 
well-chosen, so that the whole back-and-forth procedure 
induces a contracting error dynamic, one can ensure the 
convergence of the estimator towards the initial state. 

In particular, Ramdani-et-al [4], [5] have considered (for 
the case of wave equations) the theoretical study of this 
problem applying techniques borrowed from semigroups 
theory. 

Here, we consider a similar problem to [5] for a wave 
equation, where the initial conditions are known (are zero) 
and we are rather interested in identifying an unknown source 
term q{x). Let also T > 0, a; e M and let g € i?^(0, 1) n 
Hq{Q, 1). We consider the following system 

utt - u^x = qix) coa{u}t), (t, x) e (0, T) x (0, 1), 

u{t,0) = u(i,l) = 0, te (0,T), 

u{0,x) =uti0,x) = 0, x€ (0,1), 

y{t) = u,it,0),t€{0,T), 

(1) 
where {u, ut) represents the state of the system and y is the 
output. The term cos{ijjt)q{x), where a; is a fixed (known) 
frequency and q{x) is unknown, is considered to be an 
external force which varies harmonically. 
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For any q E iJ^(0, 1) n Hq{0,1) there exists a unique 
solution u e Ci([0,T];i^o(0, 1)) nC2([0,T];L2(0, 1)) to 
^ and u^{t,0) e iJi(0,T) (see [6, Remark 2]). It is 
moreover well-known from the work [6] of Yamamoto that 
this problem is well-posed in the sense that one can retrieve 
the source term q{x) from the measurement of y on the 
time interval (0, T) if T is large enough. Our aim here is 
to propose well-chosen observers which allow, using back 
and forth procedure, to retrieve q from y in the minimal 
observation time. More precisely, we prove the following 
result 

Claim 1 We can construct efficient observers for which 
the back-and-forth algorithm is convergent and which allow, 
using the measurement output y{t) over the time-interval 
(0,2), to reconstruct the unknown source term q{x). 

Remark 1: Let us point out that since the spatial domain 
is given by (0, 1), the minimum observability time is given 
by T = 2. 

Note that, whenever the whole initial state 
(m(0, .), ut(0, .), q) is unknown, system ([T]) is not observable. 
In order to realize this, one can consider the simpler case 
where the source term q is given by only the two first modes 
of the wave equation, q{x) — qi sin(7ra;) + q2 sin(27r2;), 
where qi and q2 are the unknown scalars. In this case, 
([T]i becomes equivalent to two independent oscillators with 
different frequencies and with two unknown source terms 
qi cos(a;i) and q2 cos(ajt). Moreover, the output is given by 
a linear combination of the position of the oscillators. This 
six dimensional system with one output is not observable. 
However, if we know the initial state of the oscillators, the 
two parameters qi and (72 become identifiable. 
The back-and-forth estimator allows us to take into account 
this knowledge of the initial state (w(0, .),Uf(0, .)) = (0,0) 
of the wave equation. On the contrary, if we had only used 
a forward observer, we would have lost the information on 
the initial state of the system and therefore there would 
have been no reason for the observer to converge to the real 
parameters. 

We prove the convergence of the algorithm using Lyapunov 
techniques and LaSalle's principle. One of the main 
difficulties comes from the fact that the precompactness 
of the trajectories is not ensured since we deal with an 
infinite dimensional system and we also have to use sharp 
mathematical estimates. 

The paper is organized as follows. 
In Section 2, we prove the equivalence between ([TJ and 



another system which consists in a system composed of a 
wave equation without any source term and an oscillator 
in cascade. Here the unknown term to retrieve is the initial 
condition. Section 3 is devoted to the proof of Claim 1 and 
is divided in two parts. In the first one we study the well- 
posedness of the considered systems. In the second one we 
prove the convergence of the proposed algorithm. 

II. An equivalent estimation problem 
We begin with introducing the following system, 

Wu-w.,,=Q, (t,a;)G(0,T)x(0,l), 

w(i,0) == w(t, 1) =0, te (0,r), (2) 

w{{), x) — q{x), Wt{0, x) = 0, x e (0, 1), 



ii(i) = z2(t),ie(o,r), 

Z2{t) = -uj^zi{t) + w^{t,Q), t £ (o,r), 

zi(0) = 2/(0), Z2(0) = 2/(0), 

r(i) = 2i(i),te(o,T), 



(3) 



where (w, wj, zi, Z2) represents the state of the system and 
Y is the output. System (|2]i-([3]l is nothing but a homogeneous 
wave equation coupled with an oscillator. The unknown q is 
now the initial datum of the system. We have the following 
results, whose proofs can be found in [2] 

Proposition 1: There exists a unique solution 
{w, Wt,zi, Z2) to (|2]i-([3]l with the following regularity 

«;eC([o,r]; 1/2(0, i)nffi(o,i)), 

u;teC([0,T];i/i(0,l)), 
(zi,Z2)ei72(0,T)xiji(0,T). 
Moreover, w satisfies the following hidden property 



w^it,0)eL^{0,T). 



(4) 



and 



Proposition 2: For any T > 0, for any q E H^{0, 1), 

y = Y in H^{0,T), 

where y denotes the output of system ([T]) and Y denotes the 
output of system (|2]). 

As a consequence of Proposition l2j we will now focus 
on system (|2])-(|3]l and provide well-chosen observers for this 
system in order to prove Claim 1. 

III. Observer design 

As we have already mentioned it in the introduction, we 
are going to prove Claim 1 using back-and-forth observers. 
We first rewrite system (|2]i-(|3]l in the following way 

w} = w^, (t,x) e (o,r) X (0,1), 

u,2 = u;i^, (i,a;)e(0,r)x(0,l), 
w^{t,0)=w\t,l) = 0,t&{0,T), "-^ 

w^{0,x) =q{x), w^{0,x) =0, xe (0,1), 



Zi{t) = Z2{t),tG{0,T), 

Z2{t) = -c^2^i(i) +u'i(i,o), t e (o,r), 

zi(0)=y(0),Z2(0) = y(0), ^^^ 

r(t) = zi(i),te(o,T), 

where, from above results, 

u;ieC([0,T];i/2(o,i)ni/i(0,l)), 

and 

(z\z2) eij2(0,l) xH1(0,1). 

Let us recall that our aim is to propose an algorithm which 
allows to retrieve the unknown q from the measurement 
output Y on the time interval (0, T). Our idea consists in de- 
signing a well-chosen asymptotic observer, which ensures the 
decrease of a same Lyapunov function in the back-and-forth 
procedure. To this aim, we first make the previous system 
periodical. More precisely, we define {W^,W'^,Z^, Z"^) on 
((0, +00) X (0, 1))2 X (0, T)2 as the solution to the following 
periodical system 

' Wt^ ^W^, {t,x) e (2fcT, (2fc + 1)T) X (0,1), 
W? = W^,, {t, x) e (2fcr, (2k + 1)T) X (0, 1), 
W} = -W^, (t, x) e {{2k + 1)T, {2k + 2)T) X (0, 1), 
w? = -WL, (t, x) e ((2fc + i)r, {2k + 2)T) X (0, 1), 
W^{t,Qi) = W^{t,l) = Q, te (0,+oo), 

[ WHO, x) = q{x), M/2(0, x)^0,xe (0, 1), 



(7) 



Zi{t) ^Z2{t),te {2kT, {2k + 1)T), 

Z2{t) = -iu^Zi{t) + V7i(i,0), t e (2fcr, {2k + 1)T), 

Zi{t) = -Z2{t), t e {{2k + i)r, {2k + 2)T), 

Z2{t)=Lo''Zl{t)-Wl{t,0), 

te {{2k + l)T,{2k + 2)T), 
Zi(0) = 2/(0), Z2(0)= 2/(0), 
[ Y{t) = Z^{t). 



(8) 



for k £ N. Using previous results, one easily sees that 
lVieL^,(M+;ij2(0,l)niJi(0,l)), 
W^2eL-(M^;i/i(0,i)), 
and 



{z\z^)eH, 



loci 



X hU 



Indeed, the above periodic system represents the system (|5]l- 
^ on the time intervals {2kT, {2k + 1)T), and the same 
system in the time-reversed manner on the intervals {{2k + 
l)r, {2k + 2)T). In particular, one has 

J Y{t) = Y{t - 2kT), t e {2kT, {2k + 1)T), 

I Y{t) = Y{{2k + 2)T -t),te {{2k + 1)T, {2k + 2)T). 



(9) 



The asymptotic observer we propose is the following 

■ w} = w^, {t, x) G (2fcr, (2fc + i)r) X (0, 1), 

W? = WL, {t,x) e {2kT, (2k + 1)T) X (0, 1), 

Wl = -W^', (t, a;) G ((2fc + l)r, (2fe + 2)r) X (0, 1), 

W? = -VKi,, (t, x) G ((2fc + 1)T, (2fc + 2)T) X (0, 1), 

VKi(t,0) = 7i(ii(t)-y(i)) 

+ 7i72(^3(t)-/oV'(s)ds) tG(0,oo), 
^^(t, 1) = 0, tG (0,+oo), 
. W^(0,x) = W'^{Q,x) = Q,x£{Q,l), 

Zi{t) = i2(t) - l2{Zi{t) - y(t)), t G (2fcT, (2fc + 1)T), 

Z2{t) = -uj^Zi{t) + w^(t,o), iG (2fcr,(2fc + i)r), 
ii(i) = -i2(t)-72(ii(t)-y(i)), 

tG ((2fe+l)T,(2fc + 2)T), 

tG ((2fc+l)T,(2fc + 2)T), 
Z3(t) = ii(t), iG(0,oo), 

^l(O) =^2(0) =^3(0)=0. 

(10) 
Here 71 and 72, the observer gains, are strictly positive 
constants to be fixed as the design parameters. Before study- 
ing the well-posedness of system (|9]l-([T0]i, let us introduce 
the error dynamics. The error term being defined as the 
difference between the observer and the observed system, 
error equations are the following 

Wt^ ^W^^ {t,x) e (2fcT, (2fc+l)T) X (0,1), 

W? = W^,, (t, x) e (2fcr, (2fc + 1)T) X (0, 1), 

Wl = -W^, (t, x) e {{2k + 1)T, {2k + 2)T) X (0, 1), 

Wf = -Wl,. (t, x) e {{2k + l)r, {2k + 2)T) X (0, 1), 

1^1(^,0) = 7i^i(i) +7i72^3(i), t e (0,00), 

W^t,l) = 0, te (0,+cx)), 

W^{0, x) = -q{x), W^{0, x)^0,xe (0, 1), 

(11) 
Z,{t) = Z2{t) - 72^1 W, t e {2kT, {2k + 1)T), 

Z2{t) == -oj'^Zi{t) + W^{t,0), t e {2kT, {2k + 1)T), 

Zi{t) = -Z2{t)-j2Zl{t), 

te {{2k + l)T,{2k + 2)T), 

Z2{t)^u;^Z,{t)-W^{t,0), 

te ((2fc + l)r,(2/c + 2)T), 

Z3{t) = Zi{t),te (0,00), 

Zi(0) = -y(0), Z2(0) = -m, ^3(0) = 0. 



(12) 



. Well-posedness 
From now and until the end, we denote 



H^{Q, 1) :== {w e iJi(0, 1) s.t. w(l) == 0}. 
We have the following proposition (see [2] for the proof). 



Proposition 3: For any T > 0, for any {q'^,q^) € 
F2(0, 1) n FoHO, 1) X H^{0, 1) and any (^0,^2°, Cg) e M^, 
there exists a unique solution (wi, W2,Ci 1^2:^3) to the fol- 
lowing periodical Cauchy problem 

v} = w2, (t, x) e (2/cT, (2A: + l)r) x (0, 1), 

v^ = vl,, {t, x) e {2kT, {2k + 1)T) X (0, 1), 

v^ = -1-2, (i, x) e ((2fc + 1)T, (2/c + 2)r) X (0, 1), 

«2 = -«i^, (t, x) e ((2fc + l)r, (2fc + 2)T) X (0, 1), 

v\t, 1) = 0, v\t,0) - 71^1 (t) + 7i72e3(t), 

te {{2k + l)T,{2k + 2)T), 
v^O,x) ^q°{x), v^{0,x) = q^{x), xe (0,1), 



(13) 



6(i) = 6 W - 72^1(0, t € (2fcr, (2fc + 1)T), 
6(0 = -^'Ci W + vl{t, 0), t e (2fcT, {2k + 1)T), 

t£ {{2k + l)T,{2k + 2)T), 
i2{t)^u:^Ut)^vl{t,0), 

te {{2k + l)T,{2k + 2)T), 

Ut)=Ut),te{0,^), 

a(0) = e?,6(0)=6°,6(0) = e3°, 

(14) 
where k e N, with the following regularity 

v^ eC{[0,+oo);H^O,l)nH^{0,l)), 

D2eC([0,+(^);i/,i(0,l)), 

(6,6,6) e i/i([0,+oo)) X L\[0,+^)) X h2([o,+^)). 

Moreover (w^, f^,6,6,6) satisfies the following energy 
identities: for any t e (0, +00), 

l^'(i,-)li2(o,i) + b'(^,-)li2(o,i) + 7il6WP 

+ 7iw'l6WI' + 27i72^'|aii2(o,*) = 

and 

2 (l^?(*, 011,2(0,1) + k^(i, Oli^fo,!) + 7ia;*|Ci(i)|^) 

1 1 ■ A* 

+ ^7i|i^i(i,0)|^ + -7172^^^ / |6(s)l^rfs + 7ia;^|^2(t)l^ < 

'~^\\1 Ih2{0,1) + !<? ljfl{0,l) + l?ll +|S2| ), 

where C denotes a positive constant which only depends on 

'^,71,72- 

One easily sees that the well-posedness of system ( [TT| l- 
T2I1 in 



c([o, +00); ij2(o, 1) n iJi(o, 1)) X c([o, +00); i/i(o, 1)) 

X iJi([0,+cx))) X L^{[0,+oo)) X i/2([o,+oo)) 

directly follows from Proposition [3] The well-posedness of 
system (|9|l-([T0li in 

Lf^,{R+;H^{0, 1) n i?i(0, 1)) X L^,(M+; ffi(0, 1)) 
X i//,,(M+) X LL(K+) X i/L(M+) 



then follows (using also the well-posedness of system (|7|-([8]) 
in the same space). 

B. Asymptotic analysis 

The main goal of this section is to prove the following 
result, implying Claim 1 

Theorem 1: For any T > 2, 

lim W\2nT,.)^0m Hl{0,l), (15) 

n— >-+oo 

lim W^{2nT,.) = 0mL^{0,l), (16) 

lim {Zi{2nT),Z2{2nT), Z:i{2nT)) = (0, 0, 0). (17) 

Proof of Theorem [7] Let us assume that ([T5]l-([T7]i do 
not hold. Then, there exist a positive constant a and a 
subsequence ((/'(n))„>o G N^ with lim ^(n) = +oo such 

that for any n > 0, 

\{W\2^{n)T, .),W\2(p{n)T, .))|«i(o,i)xl2(o,i) + 

|(Zi(20(n)r),Z2(20(n)T),^3(2<^(n)T))| > a. (18) 

From Proposition [3] 

W^ eL°°(M+;ij2(0,l)ni/,i(0,l)), 
W'^ <EL°^{R+;Hl{Q,l)), 
(Zi,Z2,Z3) e H\R+) X L'{R+) X ij2(M+), 
and for any time t G M^, 
\Wl{t, .)li.(oa) + |W^'(^' ^^.^(o,!) + 7i|^2 WP 

Since {(W^(2(/)(n)r, .), iy2(-2^(^)y^ ,))^ „ ^ I^j. i^ 
bounded in i7^(0, 1) x Hl{Q,l), it follows from Kato- 
Rellich's theorem that there exists a subsequence of 
{<f){n))n>o, that, for convenience, we still denote (0(n))„>o, 
and there exists {W'^^^,W°°^'^) € iJi(0, 1) x L^{0, 1) such 
that 

lim W\2(l}{n)T, .) = 1^°°^^ in H^O, 1), (20) 

n— ^+oo 



(Xi, X2, X3) e L2(0, 1) X L2(0, 1) X L2(0, 1) 

be solution to 

v} =i;2, (t,a;)e (0,T) x (0,1), 
w2=i;i^, (t,x) e (0,r) X (0,1), 

i;i(<, 1) = 0, v\t, 0) = 7ixi(t) + 7i72a^3(t), i e (0, T), 
fi(t) = a;2(t)-72a;i(i),<e(0,T), 
i2(t) = -w22.^(i)_^„i(^^0), te (0,T), 
i3(t) = a:i(t),ie(0,r), 

wi(0,a;) = W^°°'i(a;), ^^(O,^;) = W°°^'^ix), xe (0,1), 
[ a;i(0) = Z^, X2(0) = ^2°°, 3:3(0) = ^3°°- 

(25) 
Finally, let us define, for any n > 0, 



(v'''\v'n ^ C{[0,T];HX0,1)) X C{[0,T];L\0,1)) 



and 



by 



{x^,x^,x^) e L^{0,1) X L^{0,1) X L'{0,1) 

«!'"(*, .):=(Tyi(20(n)r + t,.), (26) 

v^-"{t,.):^W^{2(j}{n)T + t,.) (27) 



and 



lim W%2<j){n)T, .) = W°"-^ in L^(0, 1) 



(21) 



From ([19), #i(20(n)T, 0), Zi(2(/)(n)r) and Z2{2(j}{n)T) 
are bounded, and from the fifth equation in ( [TT| i so is 
Z^{2(f){n)T). Thus, (up to a subsequence of (0(n))„>o), 
there exists (Zf^, Z|°, Zf°) e M^ such that 



lim Zi(2(?!)(n)r) = Z^"" 

n— )-+oo 

lim Z2{2(t){n)T) = Z. 

n— ^+00 



lim Z3(2c/.(n)r)) - ^3°= 

n— J-+00 



CX3 

2 ' 



(22) 
(23) 
(24) 



cc^(i):=Zi(20(n)T + t), (28) 

a;^(i) :=Z2(20(n)T + t), (29) 

x^(i)):=Z3(2,^(n)T + t) (30) 

for any t G (0,r). With such a definition, for any n > 0, 

{v^'" , v'^'" , Xi , X2 , x^) is solution to 

i;^" = «2.„^(^^^)g(0,r)x(0,l), 

t;^" = «ii", (t,a;)e(0,r)x(0,l), 

t;i'"(t,l) = 0, i;i'"(i,0) = 7ia;?W+7i72a:^W, 
ie(0,T), 

i5'(i)=4Hi)-72a;ni),ie(0,T), 

±5 (i) = -0.2^1' (i) + «!." (i, 0) , t e (0, T) , 

i^(i)=x5'(i),ie(0,T), 

wi'"(0,a;) = T4^i(2(/)(n)r,a;), x e (0,1), 
i;2^"(0, x) = T^2|20(„)r^ x), x e (0, 1), 

x^{0) ^ Z,i2^{n)T), 

x^{0) = Z2{2<j>{n)T), 

, x^m = Z3(</.(n)T). 

(31) 
From ([T9]l-(p5]l and ( (3T] i (more precisely by continuity of 
flow with respect to the initial state). 



Let now 



(z;\«2)gC'([0,r];iJi(0,l))xC([0,r];L2(o,l)), 



«i in L°°((0,T); 1/1(0,1)), (32) 

)), (33) 

lim (a:^x5,4) = (xi,X2,a;3)inL°°(0,T)3. (34) 



lim w^^" 

n— >4-oo 

1;™ ..2,n _ ^,2 • 7-oo//n T^\. r2 



n— >-+oo 



We now introduce the following Lyapunov function 

'i\W^{t,x)\^ + \W\t,x)f)dx 
+ 7ia;^lZi(i)|^+7i|Z2(i)|'),iG (0,r). (35) 

Indeed, V(t) > 0, i > and using ([TT} and ( [T2] l, one can 
compute, for any time t E (0,T), 



We recall that if W°°^^{x) = ^ Ofc cos(^^^^^J^) and 



+ 00 



+ 00 

E' 

k=l 



W°°'^{x) — ^ bk cos(^^ — ^2^^), then v^ can be expressed 



V [t, X) = 2_^ I flfe COS ( 



26fc . f (2k + l)-!vt\ ] f(2k + l)Tvx 

^ sm ^ — cos ^ — 

(2fc + l)7r V 2 yJ V 2 



di 



V{W\W\ Zi,Z2, Z3){t) = -7172^1^1 (t) I ^ < 0. 



(i,x)G(o,r)x(o,i). 



(36) 



At x = 0, we have 



The function V is positive, decreasing. Consequently there 
exists I > such that 

lim V{W\W^,Zi,Z2,Z3){t) = l. (37) 



oo 



/'{2k + l)nt 
V 2 



On the other side, from ( [26| l, ( [28) , ( |32| i and ( |35| l, one has, 
for any te (0,r), 



lim V{W\W^, W^,Zi,Z2, Z3){2(j){n)T + t) = 



{\vlit,x)\^ + \v^it,x)\^)dx 







+ 7ic^2|^^(^)|2^^^|^^(^)|2^^ (38) 

Thus it follows from ([37]) and ( [38] ) that for any time t e 
(0,T), 

V(i;\w^,a;i,a;2,a;3)(0 = ?■ 

In other words, t i-> V(u^,ti^,xi, 0:2, a;3)(t) is constant on 
(0, T) and thus 



2bk . /(2fc + l)7rt\ 1 ,„ ^, 

Multiplying this last equality by cos ( - — \^ '^ J , for any fc G N, 
and integrating on (0, 2) we obtain 

afe = 0, fc G N, 

and thus, 

W^°°'^ =0on (0, 1). (44) 

In particular, 



and 



VF°°'i(0)= 7172^3°° = 



«,^(t,0)=£6.cos((2^^il^),tG(0,r). 



(45) 



Parseval's identity then implies 

fT /.2 



dt 



V{v ,v ,xi,a;2,X3) 



7i72w2|xi(t)|2 = 0, ie(0,r). (39) 



/ |^t(i,0)|'dt> / \vl{t,Q)\' dt = 
Jo Jo 

00 

J2\b,f^\W^'^\l2^o,iy (46) 



We finally obtain 



Thus, from l|43^, we finally get 



xi EEOon (0,T). 



(40) 



and 



Consequently, from the fourth and sixth lines in ( [25| l, we 
also get 



Vy°°.i = VF°°.2 = on (0, 1) 



Zr = ^2°° = ^3°° = 0. 



a:2 = 0, X3 = ^3=" on(0,T), 



(41) 



where Z^ is a real constant. Then, using the third and fifth 
lines in ( p5] ) and, ( ]40] i and ( ]4T| we see that 

«i(t,0) = 0, i>i(i,0)= 7172^3°° on(0,T). (42) 
Consequently, (]25]l reduces to 



f ^;l=i;^(t,a;)e(0,T)x(0,l), 
^2=z;i^, (i,a:)G(0,r)x(0,l), 
t;i(t,l)=«i(t,0) = 0, iG(0,r), 

t;i(t,o) =7172^3°°, ie(o,r), 

1-1(0, x) ^W°"'^{x), v^{0,x) = W°°''^{x), xe (0,1), 
[ xi{t) = X2{t) = 0, X3{t) = Zl", t G (0,T). 

(43) 



This is a contradiction with ([TsJ and finishes the proof of the 
Theorem [T] 

IV. Numerical simulations 

In this section, we illustrate the efficiency of the above source 
estimation algorithm through numerical simulations. We consider 
the system {TJ with source term q := x — x^ , together with the 
estimation algorithm ([9)-(|10[) with initial estimate g = 0. We fix 
the observation horizon to T = 3 and consider 50 iterations of 
the estimator ([9)-(|10^. The simulations of Figures land 2 illustrate 
the performance when we have added 10% white noise on the 
measurement output and where the observer gains 71 and 72 are 
chosen to be 

71 = 1, 72 = 1/2. 

The numerical simulations have been done through a finite dif- 
ference method where the time and the space are discretized 
simultaneously. We have chosen a spatial discretization with 20 
steps (Ax = .05) and a CFL coefficient of .005 (At = 2.5e - 04). 




Iteration nimiber 



Spatial coordinate x 



Fig. 1. The source estimator q is traced after each iteration; as it can be 
seen, the estimator after 50 iterations has converged towards q. 



15 20 25 30 35 40 45 50 




Iteration number 



Fig. 2. The first plot illustrates the decrease of the L^ distance between 
q and q after each iteration and its convergence to zeros. The second plot 
illustrates the decrease of the Lyapunov function defined in {35). 
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